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Abstract 

We use results for the structure function Fl for a gluon target having nonzero 
transverse momentum square at order a s , obtained in our previous paper, to 
compare with recent HI experimental data for Fl at fixed W values and with 
collinear GRV predictions at LO and NLO approximation. 
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The longitudinal structure function (SF) Fl(x, Q 2 ) is a very sensitive QCD characteristic 
and is directly connected to the gluon content of the proton. It is equal to zero in the parton 
model with spin— 1/2 partons and has got nonzero values in the framework of perturbative 
Quantum Chromodynamics. The perturbative QCD, however, leads to a quite controversal 
results. At the leading order (LO) approximation Fl amounts to about 10 -j- 20% of the 
corresponding F 2 values at large Q 2 range and, thus, it has got quite large contributions 
at low x range. The next-to-leading order (NLO) corrections to the longitudinal coefficient 
function are large and negative at small x [l]-[5] and can lead to negative F L values at 
low x and low Q 2 values (see [5, 6]). Negative Fl values demonstrate a limitations of the 
applicability of perturbation theory and the necessity of a resummation procedure, that leads 
to coupling constant scale higher than Q 2 (see [5], [7]- [9]). 

The experimental extraction of F L data requires a rather cumbersome procedure, espe- 
cially at small values of x (see [10], for example). Recently, however, there have been pre- 
sented new precise preliminary HI data [11] on the longitudinal SF F L , which have probed 
the small-x region 10~ 5 < x < 10~ 2 . 

In this paper the standard perturbative QCD formulas and also the so called /c T -factorization 
approach [12] based on Balitsky-Fadin-Kuraev-Lipatov (BFKL) dynamics [13] (see also re- 
cent review [14] and references therein) is used for the analysis of the above data. The 
perturbative QCD approach is called hereafter as collinear approximation and applied at 
LO and NLO levels using GRV parameterizations for partion densities (see [15]). The cor- 
responding coefficient functions are taken from the papers [3, 1]. 

In the framework of the /^-factorization approach, which is of primary consideration 
in our paper, a study of the longitudinal SF F L has been done firstly in Ref. [16], where 
the small x asymptotics of Fl has been obtained analytically using the BFKL results for 
the Mellin transform of the unintegrated gluon distribution and the longitudinal Wilson 
coefficient functions for the full perturbative series has been calculated at asymptotically 
small x values. In this note we follow a more phenomeno logical approach [18] where we 
analyzed Fl data in a broader range at small x and, thus, we use parameterizations of the 
unintegrated gluon distribution function $ g (x, kj_) (see Ref. [14]). 

A similar study has been already done 2 in our paper [18] using previous HI data [21]. 
The recent HI preliminary experimental data [11] is essentially more precise, that stimulates 
the present additional study. 

1. The unintegrated gluon distribution $ g (x, k 2 ^) (f g is the (integrated) gluon distribu- 
tion in the proton multiplied by x and k± is the transverse part of the gluon 4-momentum 



is the basic dynamical quantity in the /c^-factorization approach 3 . It satisfies the BFKL 
equation [13]. 

2 Note that the studies of the Fl structure function in the framework of the /cy-factorization have been 
done also in [19, 20]. 

3 In our previous analysis [17] we have shown that the property k 2 — —k"j_ leads to the equality of the 
Bjorkcn x value in the standard renormalization-group approach and in the Sudakov one. 
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Then, in the /cy-factorization the SF F2,l{x, Q 2 ) are driven at small x primarily by gluons 
and are related in the following way to the unintegrated distribution $ g (x, kj_): 

F 2>L (x,Q 2 ) = - dkl Y, elCi L (x/z,Q\mlkl)$ 9 (z,kl), (2) 

JX Z J -7 

where ef are charge squares of active quarks. 

The functions C^ L {x, Q 2 , m 2 , k±) can be regarded as SF of the off-shell gluons with vir- 
tually kj_ (hereafter we call them hard structure functions by analogy with similar relations 
between cross-sections and hard cross-sections). They are described by the sum of the quark 
box (and crossed box) diagram contribution to the photon-gluon interaction (see, for exam- 
ple, Fig. 1 in [17] and [18]). 

2. Notice that the /^-integral in Eqs. (1) and (2) can be divergent at lower limit, at 
least for some parameterizations of & g (x, k\). To overcome the problem we change the low 
Q 2 asymptotics of the QCD coupling constant within hard structure functions. We apply 
here two models: the "freezing" procedure and Shirkov-Solovtsov analytization. 

The "freezing" of the strong coupling constant is very popular phenomenological model 
for infrared behavior of a s (Q 2 ). The "freezing" can be done in the hard way and in the soft 
way. 

In the hard case (see [22], for example), the strong coupling constant itself is modified: 
it is taken to be constant at all Q 2 values less then some Ql, i.e. a s (Q 2 ) = a s (Ql), if 
Q 2 <Q 2 . 

In the soft case (see [20], for example), the subject of the modification is the argument of 
the strong coupling constant. It contains the shift Q 2 — ► Q 2 + M 2 , where M is an additional 
scale, which strongly modifies the infrared a s properties. For massless produced quarks, 
p-meson mass m p is usually taken as the M value, i.e. M = m p . In the case of massive 
quarks with mass m*, the M = 2rrii value is usually used. Below we will use the soft version 
of "freezing" procedure. 

Shirkov and Solovtsov proposed [23] a procedure of analytization of the strong coupling 
constant a s (Q 2 ), which leads to a new strong analytical coupling constant a an (Q 2 ) having 
nonstandard infrared properties. We are not in position to discuss here theoretical aspects of 
the procedure and use only the final formulae for the analytical coupling constant a an {Q 2 ). 
They have the following form 



a an {Q 2 ) 
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in the LO approximation and 



a an (Q 2 ) 



47T 



1 

Wo 



1 A 2 _ A 2 
ln(Q 2 /A 2 ) + h ln[l + ln^/A 2 )/^] + 2 A 2 - Q 2 ~ Q 2 'J ' 

the two first terms in the a s -expansion of 



in the NLO approximation, where (3q and f3\ are the two first term; 
/3-function and b\ = Pi/Pq- The constant C\ = 0.0354 is very small. 
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The first terms in the r.h.s. of Eqs. (3) and (4) are the standard LO and NLO represen- 
tations for a s (Q 2 ). The additional terms modify its infrared properties. 

Note that numerically both infrared transformations, the "freezing" procedure and Shirkov- 
Solovtsov analytization, lead to very close results (see below Fig. 1 and also Ref. [24] and 
discussion therein). 
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Figure 1: Q 2 dependence of F L (x, Q 2 ) (at fixed W = 276 GeV). The HI preliminary, e + p and 
e~p experimental data are shown as the black points, black and white squares, respectively 
(see [11]). Theoretical curves obtained in the kr— factorization approach with the JB unin- 
tegrated gluon distribution: solid curve corresponds to "frozen" coupling constant, dashed 
curve - analytical coupling constant, dash-dotted - "frozen" argument of the unintegrated 
gluon distribution function. 

3. As it was already noted above, the purpose of the paper is to describe new prelim- 
inary HI experimental data for the longitudinal SF F L (x,Q 2 ) using our calculations of the 
hard SF C 9 2 L (x, Q 2 , m 2 , k\) given in our previous study [17] and infrared modifications of 
a s (Q 2 ), explained above. For the unintegrated gluon distribution $(x, kj_, Ql) we use the 
so called Blumlein's parametrization (JB) [25]. Note that there are also several other popu- 



3 



lar parameterizations, which give quite similar results excepting, perhaps, the contributions 
from the small /c 2 -range: k\ < 1 GeV 2 (see Ref. [14] and references therein). 

The JB form depends strongly on the Pomeron intercept value. In different models the 
Pomeron intercept has different values (see [26]). So, in our calculations we apply the HI 
parameterization [27] based on the corresponding HI data, which are in good agreement 
with perturbative QCD (see Refs. [27, 28]). 

We calculate the SF F L as the sum of two types of contributions: the charm quark one 
F[ and the light quark one F l L : 

F L = F[ + F C L (5) 

For the F l L part we use the massless limit of hard SF (see [17, 18]). We always use / = 4 
in our fits, because our results depend very weakly on the exact / value (for similar results 
see fits of experimental data in [29] and discussions therein) . The weak dependence comes 
from two basic properties. Firstly, the charm part of Fl, F£, is quite small at the considered 
Q 2 values (see Ref. [18] for the F[ study). Secondly, the strong coupling constant very 
weakly depends on / because of the corresponding relations between A values at different / 
(see [30]). 

In Fig. 1 we show the SF Fl with "frozen" and analytical coupling constants, respectively, 
as a function of Q 2 for fixed W in comparison with HI experimental data sets (see [11]). The 
results are mostly coincide with each other. They are presented as bold and dashed curves, 
which cannot be really resolved in the figure. 

The dash-dotted curve shows the results obtained with "frozen" argument of the unin- 
tegrated gluon density. The difference between the bold and dash-dotted lines is not so big, 
that demonstrates the unimportance of the infrared modifications of the density argument. 
Below we only restrict ourselves only to the modification of the argument in the strong 
coupling constant entering the hard structure function. 

Fig. 2 contains the same bold curve as Fig. 1 and shows also the collinear results for Fl 
values. We use the popular GRV parameterizations [15] at LO and NLO approximations. 
The A/r-factorization results lie between the collinear ones, that demonstrates clearly the 
particular resummation of high-order collinear contributions at small x values in the hr- 
factorization approach. 

We also see exellent agreement between the experimental data and collinear approach 
with GRV parton densities at NLO approximation. The NLO corrections are large and 
negative and decrease the F L value by an approximate factor of 2 at Q 2 < 10 GeV 2 . 

In Figs. 1 and 2, our /c T -factorization results are in good agreement with the data for 
large and small parts of the Q 2 range. We have, however, some disagreement between the 
data and theoretical predictions at Q 2 ~ 3 GeV 2 . The disagreement exists in both cases: for 
collinear QCD approach at the LO approximation and for /c^-factorization. 

Comparing these results with Fig. 4 of Lobodzinska's talk in Ref. [11] we conclude that 
the disagreement comes from the usage of the LO approximation. Unfortunately, at the 
moment in the /^-factorization approach only the LO terms are available. The calculation 
of the NLO corrections is a very complicated problem (see [31] and discussion therein). 

A rough estimation of the NLO corrections in the /^-factorization approach can be done 
in the following way. Consider first the BFKL approach. A popular resummation of the 
NLO corrections is done in [8] at some approximation. Ref. [8] demonstrates, that is the 
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Figure 2: Q 2 dependence of F L (x,Q 2 ) (at fixed W = 276 GeV). The experimental points 
are as in Fig. 1. Solid curve is the result of the k T — factorization approach with the JB 
unintegrated gluon distribution and "frozen" coupling constant, dashed curve - the GRV LO 
calculations, dash-dotted curve - the GRV NLO calculations, dotted curve - the result of the 
GRV LO calculations with /i 2 = 127Q 2 . 

basic effect of the NLO corrections, that is the strong rise of the a s argument from Q 2 to 
<5g// = K ■ Q 2 , where K = 127, i.e. K » 1, which is in agreement with [5], [7] and [9]. 

The use of the effective argument Qlff in the DGLAP approach at LO approximation 
leads to results which are very close to the ones obtained in the case of NLO approximation: 
see the dot-dashed and dotted curves in Fig. 2. Thus, we hope that the effective argument 
represents the basic effect of the NLO corrections in the framework of the /cT-factorization, 
which in some sense lies between the DGLAP and BFKL approaches as it was noted above 
already. 

The necessity of large effective arguments is also demonstrated in Fig. 3, where we show 
the /^-factorization and collinear results for nonrunning coupling constant. Its argument is 
fixed at Q 2 = Mf giving a s ~ 0.118 (see [32]), i.e. the considered argument is larger than 
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Figure 3: Q 2 dependence of F L (x,Q 2 ) (at fixed W = 276 GeV). The experimental points 
are as in Fig. 1. Solid curve is the result of the kr— factorization approach with the JB 
unintegrated gluon distribution and ji 2 = Mf , dashed curve - the GRV LO calculations at 
f = M 2 . 

the most part of the Q 2 -values of the considered experimental data. 4 

The results obtained in the /cT-factorization and collinear approaches based on Q 2 ^ 
argument are presented in Fig. 4. In comparison with the ones shown in Fig. 1, they 
are close to each other because the effective argument is essentially larger than the Q 2 
value. There is very good agreement between the experimental data and both theoretical 
approaches. 

Moreover, we also present in Fig.4 the F L results based on the -R^o^-parameterization 
for the R = <3l/&t ratio (see [33]) (because F L = F 2 R/(1 + R)), improved in [34, 35] for 
low Q 2 values and the parameterization of F 2 data used in the our previous paper [18]. The 
results are in good agreement with other theoretical predictions as well as with experimental 
data. 

4 The study is also initiated by conversation with L.Lonnblad, we thank him. 



6 




Figure 4: Q 2 dependence of Fl(x,Q 2 ) (at fixed W = 276 GeV). The experimental points 
are as in Fig. 1. Solid curve is the result of the factorization approach with the JB 
unintegrated gluon distribution and at /i 2 = 127Q 2 , dashed curve - the GRV LO calculations 
at fi 2 = Y17Q 2 , dash-dotted curve - from the i?^ or ^-parametrization. 

4. Resume. In the framework of /c^-factorization we have applied the results of the 
calculation of the perturbative parts for the structure functions Fl and F£ for a gluon 
target, having nonzero momentum square, in the process of photon-gluon fusion [17, 18] to 
the analysis of recent HI preliminary data. The perturbative QCD predictions are presented 
also at LO and NLO approximations. 

We have found very good agreement between the experimental data and collinear re- 
sults based on GRV parameterization at NLO approximation. The LO collinear and hr- 
factorization results show disagreement with the data at some Q 2 values. We argued that 
the disagreement comes from the absence of the NLO corrections in the framework of the 
/cT-factorization. We modeled these NLO corrections by choosing large effective argument 
of the strong coupling constant and argued for our choice. The effective corrections signifi- 
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cantly improve the agreement with the HI data under consideration. 
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